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ABSTRACT. We study an urn model introduced in the paper of Chen and Wei [2], where at each discrete 
time step m balls are drawn at random from the urn containing colors white and black. Balls are added 
to the urn according to the inspected colors, generalizing the well known Polya-Eggenberger urn model, 
case m = 1. We provide exact expressions for the expectation and the variance of the number of white 
balls after n draws, and determine the structure of higher moments. Moreover, we discuss extensions 
to more than two colors. Furthermore, we introduce and discuss a new urn model where the sampling 
of the m balls is carried out in a step-by-step fashion, and also introduce a generalized Friedman's urn 
model. 



1. Introduction 

Polya-Eggenberger urn models are denned as follows. At the start, time zero, the urn contains Wq 
white balls and Bq black balls. The evolution of the urn occurs in discrete time steps. At every 
step a ball is chosen at random from the urn. The color of the ball is inspected and then the ball is 
reinserted into the urn. According to the observed color of the ball, balls are added/removed due to 
the following rules. If we have chosen a white ball, we put into the urn a white balls and b black balls, 
but if we have chosen a black ball, we put into the urn c white balls and d black balls. The values 
a, b, c, d S Z are fixed integer values and the urn model is specified by the 2x2 ball replacement matrix 
M = ( a c b d ). Urn models are simple, useful mathematical tools for describing many evolutionary 
processes in diverse fields of application such as analysis of algorithms and data structures, statistics 
and genetics, see Johnson and Kotz [9], Kotz and N. Balakrishnan [10], and also Mahmoud [11]. 

One of the most fundamental urn models is original the Polya-Eggenberger urn model [3], associated 
with the ball replacement matrix M = ( J \ ). The Polya-Eggenberger urn is a balanced urn model, the 
total number of added or removed balls is constant, independently of the observed color. A parameter 
of interest is the number W n of white balls contained in the urn after n draws. Various generalizations 
of this urn model have been considered by Hill [7], Bagci and Pal [1], Pemantle [12], Gouet [5, 6], 
Schreiber [14]; we refer to Chen and Wei [2] for a brief discussion of the aforehand mentioned works. 

This work is devoted to the study of a generalization of the Polya-Eggenberger urn model, where 
several balls are drawn at each discrete time step, their colors are inspected, and they reinserted into 
the urn. The addition/removal of balls depends on the combinations of colors of the drawn balls. 
Such urn models recently received some attention in the literature, see for example Chen and Wei [2], 
Mahmoud [11], and Renlund [13]. Chen and Wei [2] introduced a particular urn model they called 
model M , where m > 1 balls are drawn from the urn at each discrete time step. Say m — £ white balls 
and £ black balls have been drawn, < £ < m, their colors are noted, and the drawn balls are returned 
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to the urn together with addition c(m — £) white balls and c£ black balls. The ball replacement matrix 
of this urn model is a rectangular matrix M, given by 



M 



mc 
(m — l)c 

c 






c 

(m — l)c 
mc 



(1) 



with parameter c G N and m > 1. The rows of the rectangular replacement matrix encode the 
sampling scheme in the obvious way, the £-th row corresponds to the case of drawing a combination 
of (m — £) white balls and I black balls, < £ < m, where c(m — £) white balls and c£ black balls are 
being added. Note that in model M the drawing of the m balls occurs without replacement, in other 
words the distribution of the number of white balls in the sample of size m follows a hypergeometric 
distribution. 

Chen and Wei studied the distribution of the number of white balls W n after n draws and could 
show the almost sure convergence of W n , suitably normalized, to a continuous distribution by using 
martingales. The aim of this note is to provide further insight into the limiting distribution of the 
number of white balls by providing exact expressions for the expectation and the variance of W n , 
from which one obtains the expectation and variance of the limit law. Moreover, we also show to 
obtain in principle exact expressions for arbitrary moments of the limit law. Note that the case m = 1 
corresponds to the well known Polya-Eggenberger urn M = which is completely understood, 

using different arguments such as counting arguments, or stochastic processes - see for example 
Janson [8]. It is known that proportion of white balls after n draws is a martingale, and has a beta 
distribution as the limit law with parameters b/c and w/c. Therefore, the case m = 1 is excluded 
from our study. 

Throughout this work we use the notations [?] and {^}, where [?J denotes the unsigned Stirling 
numbers of the first kind (also called Stirling cycle numbers), and the {^} denotes the Stirling numbers 
of the second kind, respectively (see [4]). Moreover, we denote with x- the ^-th falling factorial, 
= x(x - I) . . . (x - £ + I), £ > 0, with x^- = 1. 



2. The distributional equation 

We consider the urn model called model M of Chen and Wei [2], specified by ball replacement 
matrix (1). Let W n and B n denote the random variables counting the numbers of white balls and 
black balls after n draws, n > 0. We assume that the initial numbers Bq of black balls and Wq of 
white balls satisfy Bq > and Wo > 0, and that the total number of balls at time zero To satisfies 
To = Bo + Wq > m, in order to avoid any degenerate cases. Since the urn model is balanced, 
regardless of the inspected color combination a total of m ■ c balls are added to the urn at every 
discrete time step, the total number T n of balls after m draws is a deterministic quantity, and given by 



T n = T + nmc = W n + B n , n> 0. (2) 

We are interested in the random variable W n , counting the number of white balls contained in the urn 
after n draws, n > 0. The starting point of our considerations is the distributional equation, 
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W n ^W n ^ + J2k-c-l n (W k B m - k ), (3) 

k=0 

which says that the number of white balls after n draws can be decomposed as the number of white 
balls after n — 1 draws, plus the additional balls added when the colors of the n-th draw have been 
inspected, n > 1. Here the random variables I n (W k B m ~ k ), < k < m denote the indicators of 
drawing k white balls and m — k black balls from the urn at the n-th draw, n > 1. Let F n -\ denote 
the (T-field generated by the first n — 1 draws. By (2), we have -B n _i = T n _i — W n -\, and further 



mfTt Hxrk r>m—k\ _ i | 77 i _ \ k I \rn—k) _ \ k J \ m-k J 
F{IL n (W B ) - \\J- n -!\ - -f—T - jjr—^ , (4) 

\ m ) 



< k < m, n > 1. In order to study the moments E(W^), s > 1, of the random variable W n , we 
have to derive first a distributional equation for In order to do so, we take the s-th power in (3), 
and use the fact that indicator variables are mutually exclusive. We obtain the distributional equation, 
valid for n > 1, s > 1: 

s / \ m 

W* n ( = } WU + E l) W n-l^ E k " ■ UW k B m - k ). (5) 

1=1 w k=l 

3. Results for the moments 



Theorem 1. The expected value of the random variable W n , counting the numbers of white balls after 
n draws, is given by E(W n ) = ^(nmc + T ), and the variance Y(W n ) = E(W^) — E(W n ) 2 is 
determined via the second moment 

/n-l+Ai\ (n-l+X 2 \ / m J2 m ri-1 p , T -m \ 

where the values Ai, A2 are g/vew fry 

-± + mc + T ± i ^1 + 4mc(l + c) 

Al,2 — • 

mc 

Concerning higher moments, we obtain the following recursive characterization. 
Theorem 2. The s-th moment K(W^) is for s > 1 given by 

■ n— 1 \ n— 1 



j=0 a j,s 



where the quantities a„ jS , /3 n>s are defined as 

'(?) 



£=0 ui 
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s s / \ 1 /nr 3 ]( m ) 

^=E E ^ +1 iE / £ (-iy +M ^j^ Uj j . (7) 

i=2 £=i ^ ' j=l+\-i \ j ) 

Here denotes an (unsigned) Stirling number of the first kind, and {?} denotes a Stirling number 
of the second kind. 

The result above enables us to obtain in principle arbitrary high moments, since the moment E(W^) 
can be expressed terms of the moments E(Vt / £), with < h < n — 1, < r < s — 1. Note that vast 
simplifications occur in the case m = 1, which we did not observe for m > 2. As mentioned in the 
introduction the case m = 1 is very well understood; therefore we do not go into details. Moreover, 
we obtain the following results concerning the moments of the normalized random variable W n /n. 

Corollary 1. The limits lim„_ s>00 K(W^/n s ) exist; in particular, we obtain for the expected value 
linin-j.oo E (^"-) = W jT lc , and for the normalized second moment 

nwp tj^mi (w- 1 Woc2m v t + S Gil ( ^ 

Moreover, for arbitrary s > 1 f/ze //m/f o/ normalized s-th moment can be expressed in terms of 
an infinite sum 



w/f/i /3^ )S as defined above in Theorem 2, involving the moments of the form E(WT), w/?/i < /i < oo 
£m<f, < r < s — 1. //ere \i )S denote the roots (times by minus one) of the equation 

Remark 1. As mentioned in the introduction Chen and Wei [2] proved almost sure convergence of 
W n , namely W n = ^ = nT7 ^"r Woo- Hence, we get the expectation and the variance of 
Woo by the relations E(Woo) = and V(Woo) = lim^ Our results show 

that the distribution of Woo is not an ordinary beta law, in contrast to the case m = 1. Note that the 
moments of W n do not grow very fast, since they satisfy the trivial bounds E(W*) < (nmc + Tq) s . 
Hence, by Carleman's condition the limit law is uniquely determined by its moments E(W£ 3 ), 
which are given by E(W^) = lim^oo 

4. Determining the structure of the moments 
In order to prove Theorems 1 and 2, we need the following result. 
Lemma 1. The moments E(W^) satisfy the recurrence relation 

E(W°) = a n _i, s • E(W^_i) + Pn-i, 8 , n,s>l, (8) 
where the quantities a ns and fi na are as defined in Theorem 2. 
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Our starting point for the proof of Lemma 1 is the distributional equation for W% (5), and we take the 
conditional expectation with respect to T n -\. This leads to 



a s \ III, 

1=1 ^ ' k=l 



m (W n -i\ /T n _i-W„_i\ 

I \ k )\ m-k ) 



V m I 



where we have used (4) and (2). In order to simplify the stated expression we have to use some 
combinatorial identities. We convert the ordinary powers of k into falling factorials using the Stirling 
numbers of the second kind, 



k=l ^ 



where x- = x{x — l)(x — 2) . . . (x — (k — 1)), rn > 1. We have 



k=l k=l j=l ' j=l k=j 

and consequently obtain 



m— j 

V k A m-k J 1 7 J n_1 ^ V fe I \ m - j - k 

k=l v 7 v 7 j = l w ^ fc=0 V / v / 

Next we use the Vandermonde convolution formula 



E 

fc=0 

in order to obtain the expression 



r \ I s \ i r + s 
k J \n — kj \ n 



This leads to 



E*' . ' ^ » j = E VLfe 1 = E 



fe=l V m ) j=l 



Next we convert the falling factorials into ordinary powers, and obtain 

t-iy-'wj-i. 



wi.. = E 

i=l 

Hence, we have 



3 

3 
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5 CH- 1 (V)- g {;} (Y) I? 



■re 



V 3 ) 3=1 KJ J \ 3 > 1=1 

j+i-e-i (jKj) 



E<-^ E (-*) j 



=1 j=e+i-i \ j i 



Note that the result above is an, explicit expression for the ^-th moment of a hypergeometric distributed 
random variable with parameter W n -i, T n _i, and m. This implies that 



£=1 ^ ' k=l \ m ) 



^-i+eCVe^^ e (-d^- 1 ^^-f' • 

£=1 ^ ' i=l j=£+l-i V j / 

The stated result now easily follows by taking the expectation on both sides. 

Remark 2. Note that in the case c = 1 a simpler expression exists for the factorial moments E(Wn) 
K(W n (W n — 1) . . . (W n — s + 1)) of W n , and consequently also for the ordinary moments of W n . 

s fs\ (m\ s s /s\ (m\ (s—£\ t£\ 

E(wi) = E (wti) -eW+e ^(w s ^) e [e)[ i:.\ )U - 

1=0 U j i=l t=i \ £ ) 

Next we use Lemma 1 in order to prove Theorem 2; we have 

E(W£) = a n _ M • ECW*.!) + /3„_ M , n > 1, a > 1. 

Let e n s be defined as 



We have 



en, 8 = n li , n>0, with e 0jS = E(iy s ) = ^ S . 
llj=o a J> 



, Pn—l,s ^ -, 

llj=o a J> 



and consequently, 



n— 1 ^ n— 1 



1 ljf=0 "j> ^=0 1 lj=0 a J> £=0 1 lj=o a j,s 

which implies the stated result. 
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In order to obtain the result for the expected value and the variance (second moment), as stated in 
Theorem 1, we proceed as follows. In the case of the expected value we observe that fyi = 0, and 
that 

cm Tj + cm Tj+% 
a jA = 1 + TfT = — Tf, = > 

since Tj + \ = Tj + mc, the total number of balls contained in the urn increases by mc after each draw. 
Consequently, 

n-1 n-1 

n = n 3+1 - 



3=0 j=0 



Tj To 



and the stated result follows. 



Remark 3. As already mentioned before W„ = W n /T n Woo- This can easily be seen as 
follows: we readily note that for s = 1 it holds 

E(W n \T n -i) = W n - 17 ^-, n > 1. 

J-n-l 

Thus, W n = W n /T n is a positive martingale with respect to the filtration F n , as previously observed 
in [2], which directly leads to the proof of the almost sure convergence of W n . 

In order to obtain the variance V(W n ) = E(W^) — ~E(W n ) 2 , we study the second moment E(W%), 
given by 

• n— 1 \ n— 1 

P£,2 



(71,-1 \ ft — 1 

j=0 J 1=0 llj 



=0 a i,2 



The value a„ 2 is given by 



- 1 ■ 2cm ■ c2 ^) - Tn+ 2 (r n -l) + c 2 m(m-l) 
a ™' 2 " + T n + rr») " T„(T B -1) 

(n + 2 + ^-)(n + + 1 - J- 

v cm £ v cm £ m 

(n+ ^)(n + 

V cm / \ cm / 

We factor the numerator of a n ^ by determining the zeros of the quadratic equation with respect to n 
and get 

_ (ra + Ai)(ra + A 2 ) 
an > 2 ~ ( n + HL) (n + 2W)' 

\ cm ' \ cm / 

with Ai, A 2 as stated in Theorem 1. Concerning f3 n ^ we have 

a tww i2f m ( 2) A Woc 2 m / m - 1 \ W c 2 m T n - m 

= nw n ) C (_ - = — (1 - _j = — • ^- T . 

This readily leads to the stated exact result for the second moment. 
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4.1. Asymptotic expansions. We finally turn to the proof of Corollary 1. In order to get the results 

I" 



for the limit linin^oo E(W 2 )/n 2 we proceed as follows: first we write n?= o a i>2 in terms of Gamma- 



functions, 

B ff Q mrr) _ r(n + A 1 )r(n + A 2 )ra)r(^; 

Using Stirling's formula for the Gamma function 



and the fact that Ai + A2 = 2 + 2t ^°. 1 , we obtain the asymptotic expansion 

«-i r / To ip/ To-i N 1 



The asymptotic expansion of n?=o a i: 2 also implies that the sum 5Z"=o — * s convergent for 



rn— 1 ,. „|..„ ... »u, „ — l 

n — >• 00 by the comparison test, 

n— 1 ^ re— 1 .. 

for some suitable constant K > 0. 

In order to provide the asymptotics of E(W*) we proceed by induction. Before we actually do so, we 
derive an asymptotic expansion of Ilj=o a j>* We nave 



(m\ sps e(m\(Tj-i\ v-s Jlfm^fTj-ts 



SI- 



SI- 



v^s J/mWjmc+T -ft 
2^=0 6 Uil s-^ / 



'(mc)-IB = i(j + a S^) 
Let A^ s denote the roots (times minus one) of the equation 



E t < I m\ fx ■ mc + T - t \ t-t 



(mc) 
We get 



' (mc)- TlUti + ^) nuu + ^ 

and consequently 



ff„ -n £(!L±^fSgi 
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By Stirling's formula we obtain the asymptotic expansion 

J — u 77, mc t — 1 

Let [rr fc ] denote the extraction of coefficients operator. Since 



m\ ( x ■ mc + To 
s-l 



(mc) s V s! (s — 1)! / mc 



we obtain 



n— -n4<Sr( i+ < 

jf=0 £=1 v ' ; 

Now we assume that E(W*) = k s ■ n s (l + O (^) , being true for s = 1, 2. By Theorem 2 we have 

n— 1 \ n— 1 r> 



n 



<=1 V ^=0 A -1.7=0 



By our induction hypothesis, E(W^) = k& • £ k (l + C(|) for sufficiently large I, 1 < fe < s — 1, 
and consequently s satisfies /3^ s < • ^ s_2 , where K denotes a sufficiently large constant, only 
depending on s. Using our result for n*Lo a J>' tne sum S"=o tJ— 21 — i s convergent according to 

the comparison test with ji- Hence, E(W^) = k s ■ n s (l + O(^)), and we have proven our 

stated results. 

5. The case of three or more colors 

The urn model M can be readily generalized to r > 2 different colors. As in the case of two colors m 
balls are drawn at random from the urn, say ki balls of color i, 1 < i < r, their colors are noted, and 
they are returned to the urn together with c • k{ balls of color i. Let X„ = (X ni i, . . . , X n>r ) denote the 
random vector counting the number of balls of type i contained in the urn after n draws, with initial 
values Xo = (Xq^, . . . , Xq >t ). The urn is again balanced, so the number of balls after n draws is 
given by T n = To + nmc, with To = Yn=i Xq,%- One gets the distributional equation 

(d) 



X n l = , X„_ 1 + Y, c-k-I n ([[X^), 



felH \-k r =m l=\ 

k e >0 
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where k = (k\, . . . , k r ), and I n (rii=i ^e l ) denote the indicators of drawing kg balls of colour I, 

1 < £ < r at the n-th draw. Note that each individual random variable X n ^ satisfies X n ^ = 
W n , where W n denotes the previously considered random variable from the two color case. The 
distributional equation above can be used to study the mixed moments of X n ,i similar to the results 
of Theorems 1, 2, and Corollary 1. We refrain from going into details since the resulting expressions 
get very involved; we only mention our findings for the covariance of two different colors with 
r > 3 and 1 < i < j < m. 

r™,( Y Y \ — \ n J\ n J Y Y { * ^ rnc> Y Y 

*-OV^A n .j, -A n ,j) — j-: r n _i -A-OA^Oj 7^0 ^0,i A 0.j- 

In— 1+— U-\ In— 1+-" — \ i„ 

( mc \ [ mc ) U 

V n ) \ n I 

Moreover, in the limit we obtain 

T \tvTq-1- 



r X nji X nJ ._ v v ( T{^)T{^) {mc) 



lim Cov( '-, ^-) = X OA Xoj i 7PT 

n^co n n ' ' J \ r(Ai)r(A2) Tq 

with Ai, A2 as given in Theorem 1. 

6. Urn model R- drawing with replacement 

We consider another urn model which we call model R. It can be considered as a variant of model M. 
The dynamics of model R are identical to model M: m > 1 balls are drawn from the urn containing 
balls of colors black and white. Their colors are inspected, and they are returned to the urn. According 
to the observed colors we add new balls: if I black balls and m — I white balls have been observed, we 
add c(m — £) white balls and I black balls; the ball replacement matrix coincides with (1). The main 
difference in model R is that the sampling of the m balls occurs with replacement, i.e. the m balls 
are drawn one by one from the urn, the colors are observed, and then put back into the urn. Hence, 
the distribution of the number of white balls in the sample of size m follows a binomial distribution 
instead of a hypergeometric distribution in model M . The distributional equation for the number of 
white balls W n after n draws is identical to (3), 



k=0 

but the distribution of the indicator variables I n (W k B m ~ k ) changes to 



, k m _ k . . , iX) W n-l B n-l {Tl W n-l( T n-l ~ W w _i) m k 

,[W B ) = l\J- n -l\ = — = — , (10) 

± n-l 1 n-l 

< k < m, n > 1. 

Concerning a limit law for the number of white balls W n after n draws, we can extend the results for 
model M of Chen and Wei [2] to model R. 

Theorem 3. The random variable W n = W n /T n is a positive martingale with respect to the natural 
filtration T n , "E,(yV n \F ri -i) = W n _i. Consequently, 

Furthermore, for fixed Wq, Bq, c and m, Woo is absolutely continuous. 
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Moreover, the moments of W n satisfy recurrence relations similar to model M. 

Theorem 4. The expected value is given by E(W n ) = ^(nmc + Ifa) and the variance Y(W n ) 
E(W^) — E,(W n ) 2 in terms of second moment, 

where the values p,\, ^ ore given by m t 2 = To+m ^ c ^ _ The s . m moment satisfies the recurrence 
relation 

■ n— 1 \ n— 1 r 

0* . 



(ft— 1 \ /£— _L 

n^i-Ko+E 



where the quantities J ns , $n,s are defined as 

fcO n 3=2 «=j V 7 J ™ 

In the following we present recurrence relations for the limits of the moments of the normalized 
random variable W n /n, with lim^oo E(W*/n s ) = (mc) s E(W^ ). 

Corollary 2. The limits of the normalized moments K(W^/n s ) exist, , and satisfy 

with 7j s , S^ s as defined above. Here the \x^ s denote the roots (times minus one) of the equation 

(^E(:)(T)^(x. m c + T r« = n ( x + «„). 

v > e=o v 7 v 7 i=i 

In the following we first sketch the proofs of Theorem 4 and Corollary 2. Since the proofs are very 
similar to the proofs for model M, we will be very brief. Then we discuss the proof of Theorem 3. 

6.1. The structure of the moments. Our starting point is again the distributional equation for W n , 
which leads to a distributional equation for W*. We take the conditional expectation with respect to 
F n -i, and obtain 



e=i ^ 7 fc=i ^ 7 



1 n-l 



The sum appearing on the right hand side is the ^-th moment of a binomial distributed random variable 
with parameters m and W n -x/T n ^\. We get the result 



m-- 



V k ) T m , ~ ^ I 7 I "" t j 

fc=l V 7 3=1 UJ 1 n-l 
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This implies that 

E(W*|.F n _l) = Tn-l.sWn-l + S n - hs , 

with 7„ ;S , 5 ntS as stated in Theorem 4. This recurrence relation can be solved similar to the proof of 
Theorem 2 for model M, and the stated result for E(W^) follows. Moreover, we easily obtain the 
stated formula for the expected value, and also the result for the second moment of W n . The results 
for the higher moments can be obtained similar to model M. Concerning the asymptotic expansions 
one can proceed similar to the proof of Corollary 1 ; we omit the details. 

6.2. Martingales and absolute continuouity. Since 5 n -i 5 i = we obtain for W n = W n /T n 

E(W„| Fn-l) = 

Hence, W n is a martingale. Since it is a positive martingale, it converges almost surely to a limit Woo- 
Note that, one has W n ~ J^; provided W n — > oo, and thus we can obtain the moments of Woo via 
the moments of 

Concerning the absolute continuouity of the distribution of Woo one can adapt the argumentation 
of Chen and Wei [2]. For the convenience of the reader we outline the main steps, quote the main 
results of [2], and present the new ingredient of the proof for model R in Lemma 2. Let (Q, P) be 
the probability space. In order to prove that the distribution of Woo has a density, one introduces a 
sequence of events {VLe)^>\ such that C Qe+i, P-fU^f^} = 1. Then, Woo is restricted to and 
it is shown that it has a density fy. The proof is then finished by proving that / = lim^oo fe exists 
and that / is the density of Woo- 

Proposition 1 (Chen and Wei [2]). Let^^^ibe a sequence of increasing events such that PjU^r^} = 
1. If there exist nonnegative Borel measurable functions (fe)e>i such that P(f^ H W^ 1 (-£>)) = 
in fl{x)dx for all Borel sets B, thenf = lim^oo fl exists almost everywhere, and f is the density of 

Woo. 

In order to construct the sequence of events (Q^^i one can follow [2]: 
Proposition 2. For fixed Wq, Bq, c and m let 

= {uj : We(u) > cm and Bg(oj) > cm}, £ > 1. 
Then, £l e C and¥{jf =1 9,i} = 1. 

In order to show that Woo has a density by restricting Woo to Q^, it suffices to show that the restriction 
of Woo to j = {uj : We(ui) = j} has a density for each j , with cm < j < T^\. For this, the 
following result is needed, which is main new ingredient in our argumentation (compared to model 
M). 

Lemma 2. The sum YliLo P{Wn+i = j + ck\W n = j + c(k — i)} satisfies 

in 

nWn+l = j + 0k\W n = j + C(k - i)} 
1 m m—i / \ / -\ 

= ^E (7) r 7 ') c* + <h - m-i - <k - i)r-^. 
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Consequently, for a fixed I it holds for all n > £, cm < j < and k < m(n + 1), and a suitably 

choosen constant k > 0. 

y^nwn+i = j + ck\w n = j + c (k - *)} < i - - + 4- 

i=0 

Remark 4. Note that the corresponding result of Chen and Wei for model M (Lemma 4. 1 and Lemma 
4.2 in [2]) can be extended and largely simplified noting that 



E 

i=0 



j + c(k - i)\ (T n - j - c(k - i) 



m — i 



m rn—f 

E^'E 

/=0 i=0 



j + c(/s - f 
i 



^ [/](-!) 



-j7-j- c ( fc - i ) N \ 
V m-i-t ) 



n 



Proof of Lemma 2. One has 



j[>{W n+ i = J + cfc|W n = j + c(fc - i)} 

i=0 

= ^ E (7J (i + c{k " i))i(Tn " j " c{k " <))r 



j=0 



£ (7 )C* + 



£=0 



m — i 



T^-j-c(k-i)) 



m—i—l 



r Pr, 
' n 



m m—l 

E^E 

£=0 i=0 



m \ m — i 



{j + c(k-i))\-j-c(k-i)) 



m—i—i 



such that 



Using T n = nrac + To we get further 

(m - + cfc) + c _ (m - + cfc) + c + q/J_\ 
T n nmc n 2 

Since (m-i)U+ck)+c > for < /c < m(n+l), W + cm < j < the stated result follows. □ 

The proof of Theorem 3 can be finished by combining the results of Propositions 1, 2 and Lemma 2; 
it is identical to the proof in [2] (Proof of Theorem 4.2) and therefore the details are omitted. 
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7. Outlook 

An interesting line of research would be to extend the results for models M and R to non-balanced 
urns with ball replacement matrix given by 



ma 
(m — l)a 

a 





b 

(m — 1)6 
mb 



(11) 



with a, b € N. In contrast to the balanced versions it seems much more difficult to obtain closed form 
expressions for the moments of W n . 

One can analyze a generalized Friedman's urn, with ball replacement matrix 



M 





c 

(m — l)c 
mc 



mc 
(m — l)c 

c 




(12) 



and parameters c G N, and m > 1. If, say, m — I white balls and i black balls have been drawn from 
the urn, < I < m, then the drawn balls are returned to the urn together with addition c£ white balls 
and c(m — I) black balls. It is possible to look at drawing the m balls without replacement, which 
we call model FM, or to look at drawing the m balls with replacement, or simply model FR. In any 
case, the distributional equation for the number of white balls W n after n draws is given by 



(d) 



V = J W n . x + £(m - k) ■ c ■ UW k B m ~ k ), 



k=0 



with 



( (W n -i\ (T n -!-W n -i\ 
\ k A m-k J 

V m I 

1 n-l 



»{I n (W K B 



k T}m—k\ 



l|J" n _l} = < 



m—k 



model FM 



model FR 



< k < m, n > 1. We can obtain the moments of W n for both models FM and FR, similar to our 
previous results for models M and R. In particular, the expectation and the variance can be obtained; 
we get for example the result 



(mc ) 2 (™) + mcT n + (T - mc)W 



mc(n — 1) + Tq 



mc(n + 1) 



, mc / T 



mc = Tq, 



n > 0, being valid for both models FM and FR. In contrast to the generalized Polya urn models 
M and R the simple martingale structure of W n /T n is not present anymore. However, one can find 
values (p n and i[> n such that M n = if n W n + ip n is a martingale, 
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¥.{M n \F n -i) = M n -1, 

with ip n and ij) n being given by 




It seems to be possible to derive a limit law using the martingale central limit theorem, or to apply the 
methods of moments in order to obtain moment convergence. The study of the generalized Friedman 
urn models FM and FR will be the subject of a forthcoming companion work. 
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